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COMPUTER PROGRAM DOCUMENTATION FOR THE DYNAMIC ANALYSIS 


OF A NONCONTACTING MECHANICAL FACE SEAL 

by 

B. M. Auer and I. Etsion 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland , Ohio 44135 


SUMMARY 

A computer program is presented which gives a numerical solution to a set 
of second order ordinary differential equations which describe the dynamic be- 
havior of a noncontacting mechanical face seal. The equations were based on an 
analytical solution of the Reynolds equation for the fluid film pressure com- 
bined with a numerical solution, using time integration techniques, of the 
equations of motion of the flexibly mounted seal ring. 

The face seal consists of two discs, a seal seat attached to a shaft and 
rotating with it and a seal ring attached to a stationary mount. The seal ring 
is flexibly supported to allow self alignment and tracking of the mating seal 
seat. Program input includes variables describing the geometry and the physical 
characteristics of the face seal and seal operating conditions. Output from the 
program includes velocities and displacements of the seal ring along and about 
the axes of an inertial reference system. 

Five appendixes are included, containing the nomenclature (A), the diction- 
ary of Fortran symbols (B) , the computer program (C), the program flow chart (D) , 
and a sample problem (E) . 


INTRODUCTION 

A dynamic analysis of a noncontacting, coned, mechanical face seal is pre- 
sented in Ref. [l]. The seal consists of two discs, a rotor which is attached 
to the shaft and rotates with it (see Fig. 1) and a stator which is flexibly 
supported to allow self alignment and tracking of the rotor.. The stator, or 
seal ring, can move axially and tilt about two orthogonal diameters while its 
circumferential rotation is prevented by anti-rotation locks. The motion of 
the flexibly mounted stator is controlled by the forces and moments acting upon 
it, by its mass inertia, and by the stiffness and damping properties of the 
seal system. 

The dynamic analysis in Ref. [l] is based on a numerical solution of the 
equations of motion of the stator. For this purpose, a computer program was 
prepared. Input data for the program describes the geometry and the physical 
characteristics of the face seal, as well as its operating conditions. These 
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are the inner and outer seal ring radii, the seal ring mass and radius of gyra- 
tion, the coning angle, the external moment acting on the seal ring, the spring 
constant of the flexible support, the radius of the springs location, the seal 
design clearance, the inside and outside pressures, and the shaft speed. The 
computed results at any time t, include the seal ring axial and angular dis- 
placement and the minimum clearance between the ring and the seat. The first 
time derivatives of the ring displacements are also computed. 

The object of this report is to describe the equations of motion of the 
face seal and the method of solution of these equations, to describe the input 
and output of the program, and to document the computer program. 


STATEMENT OF THE PROBLEM 


The program solves the set of three second order ordinary differential 
equations, (1), (2), and (3) for a x , ay, and z as a function of time. The 
full set of equations, (1) through (28), is listed below, and this is followed 
by a brief description of their derivation. 


«X = (Mx + Me)/ 1 

a) 

«y = (My/I 

(2) 

= F z /m 

(3) 

I = \ mR 2 
2 g 

(4) 

1. + Ri 

Rm = 2 

(5) 

M x = Micos(ij>) - M2sin ^ 

(6) 

My = M^sin(i{/) + M 2 COS ^ 

(7) 


2tt „ 1 




PR 2 cos 0 dRd6 - f R 2 y — 

r 0 •'Ri r ° 


y r 2tt /• i 

/ 

o J Ri 


PR sin 0 dR d 0 


( 8 ) 


(9) 
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PR dR d0 - F 


( 10 ) 


C n I" 0(1 - Rj ) "1 

F c ■ K f + Pi - TT J U. - R t ) <P ° ' Pl) J 


( 11 ) 


P = P d + P s If (P < 0), then P = 0 


( 12 ) 


P d = 


(1 + Z) 3 


(H 


“^) YR m sin 6 " z " Y ^m cos 6 




Hf(H„ + H) 

P S = p o " (P 0 - p i> A 


(13) 


(14) 


Y = Va 2 + a 2 (15) 

"x y 

Y = (a x a x + “y a y)/Y (16) 

'f' = tan _1 (a y /a x ) (17) 

i> = (a y a x - a x a y )/y 2 (18) 

e = Y/(l + Z) (19) 

6 = 0/(1 + Z) (20) 

H = 1 + cR cos(0) + 6(R - R i ) (21) 

Hi = 1 + eRiCOs(0) (22) 

H Q = 1 + e cos 0 + 6(1 - Ri) (23) 

H m = (Hi + H 0 )/2 (24) 

A = (1 - R)/jH m H 2 (l - Ri)] (25) 
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H MIN1 

= 1-0% 

(26) 

H MIN2 

= 1 - e + 6(1 - R i ) 

(27) 

h min 

= min(h min1 , h min2 ) 

(28) 


The theoretical model is shown in Fig. 2. The seal seat is parallel to 
the xy plane of an inertial reference system xyz and is rotating at a 
constant angular velocity a) about the z axis. The primary seal ring has 
three degrees of freedom. It can move axially along the z axis and tilt 
about the x and y axes. A rotating coordinate system 123 coincides with 
the principal axes of the ring so that axis 3 is perpendicular to the plane 
of the ring. The coordinate system 123 rotates in the inertial reference xyz 
so that axis 1 always remains in the plane xy and axis 2 is directed to the 
instantaneous point of maximum film thickness* Thus, the orientation of coor- 
dinate system 123 in the inertial reference xyz can be defined by the nuta- 
tion angle y* measured from axis z to axis 3, and the precession angle ^ 
measured from axis x to axis 1. 


Equations (1), (2), and (3) are the general equations of motion of the 
seal ring where a x and oty are^ the accelerations of the ring about the x 
and y axes, respectively, and Z is the axial displacement acceleration. 
Equations (6) and (7) transform the ring moments from the rotating 123 coor- 
dinate system to the xyz inertial system. Equations (8) and (9) are used 

K „2 .. C 0 


to calculate the moments Mt and M 0 where 


1 R sp Y 




is the restoring 


moment provided by the flexible support. Equation (10) calculates the axial 
force, F z , and Eq. (11) calculates the closing force, F c , contributed by both 
the flexible support and the pressure balance of the seal. Equations (12), (13), 
and (14) calculate the fluid film pressure distribution, where P s is the hydro- 
static pressure component and is the hydrodynamic and squeeze pressure 

components. Equations (21), (22), (23), and (24) give the film thickness dis- 
tributions H, H^, H 0 » and and Eqs • (19) and (20) give the tilt and the 

coning parameters e and 6. 


The velocities a x , ay, and Z and the new displacements a x , ay, and Z 
are found from the solution of the system (1), (2), and (3) to give the trans- 
ient dynamics of the seal. From these values, the nutation y and the pre- 
cession ^ and their derivatives can be found by using Eqs. (15), (16), (17), 
and (18). After each time step, the minimum film thickness, is calculated 

from Eqs. (26), (27), and (28), and is tested against a failure criterion. 

Since contact can take place on either the inner or the outer radius of the 
ring, the minimum film thickness is set equal to the smaller of or 

h MIN2* 


METHOD OF SOLUTION 

The program solves the set of three second order differential equations, 
(1), (2), and (3), where initial conditions are given for a x , ay, Z, a x , a y . 
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and Z. These equations were rewritten in terms of the six new variables 
Y 2 > Y 3 , Y 4 , Y 5 , and Y 6 as the following set of six first order equations, 
with initial conditions given for Y-^, Y 3 , Y 3 , Y 4 , Y 5 , and Y^. 


*1 - 

. 

«x 

(29) 

y 2 - 

a y 

(30) 

y 3 = 

« 

Z 

(31) 

*4 = 

(M* + M e ) /I 

(32) 

y 5 = 

M y /I 

(33) 

y 6 - 

F z /m 

(34) 


The dependent variables are given by Eqs. (35) through (40) as follows: 


Y 1 = 

a 

X 

(35) 

Y 2 = 

a y 

(36) 

Y 3 = 

z 

(37) 

y 4 = 

a x 

(38) 

Y 5 = 

a 'y 

(39) 

y 6 = 

z 

(40) 


To find the solution of Eqs. (29) through (34), the program uses the GEAR 
differential equation solver software package in the Lewis Research Center’s 
International Mathematical and Statistical (IMSL) software library of scientific 
Univac 1100 Fortran V programs. The Lewis Research Center subscribes to this 
software package which is supplied by the International Mathematical and Sta- 
tistical Libraries, Inc. software company in Houston, Texas. 

The GEAR package finds the numerical solution of initial value problems 
for systems of first order ordinary differential equations. This package has 
the option of using two basic methods of solution, either an implicit multistep 
variable order Adams predictor corrector method or an implicit multistep back- 
ward differentiation formula method, also known as Gear’s stiff method jjf]. 

Let Y x , Y 2 , Y 3 , Y 4 , Y s , and. Y 6 be denoted by Y and let Y-^ Y 2 , Y 3 , Y 4 , 

Y 5 , and Yg be denoted by Y. Then the call to the GEAR package is CALL DGEAR 
(N 6 , FCN , FCNJ, T, DELH, Y, TEND, TOL, METH, MITER, INDEX, IWK, WK, IER) , where 
the parameters passed in the call statement are defined as follows: 
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N6 - Input giving the number of first order differential equations to 
be solved. 

FCN - User supplied subroutine for evaluating the function Y. 

FCNJ - User supplied subroutine to compute the N6xN6 Jacobian Matrix of 
the partial derivatives of Y with respect to Y. FCNJ is needed 
only when miter =1. 

T - The independent variable. T is input on the first call only, and 

on the first call passes the initial value to the program. T is 

output on subsequent calls and passes to the user the current value 
of the independent variable at which integration has been completed. 

DELH - On the first call only, DELH is input and contains the next step 

size in T to be used. On subsequent calls, DELH is output and con- 
tains the step size actually used. 

Y - The__dependent variable, Y. On input, Y supplies the initial values 
of Y. On output, Y contains the value of Y computed at TEND. 

TEND - Input giving the value of T at which a solution for Y is desired. 

TOL - Input giving the relative error bound. 

METH - Input indicating whether the Adams Method or the Backward Differen- 
tiation Formulae Method is to be used. 

MITER - An input indicating whether the chord method or functional itera- 
tion is to be used and whether a finite difference or analytic 
Jacobian is to be used with the chord method. 

INDEX - Both an input and an output parameter. Index tells whether this 
is the first call to DGEAR and whether the integration was 
successful. 

IWK - An integer work vector. 

WK - A real work vector. 

IER - An output error parameter telling the user whether the integration 
was halted due to an error. The value of IER indicates the error 
encountered, as described on the following pages. 

The GEAR package is called from the MAIN program and before the first call, 
the parameters N6, DELH, TOL, METH, MITER, and INDEX are assigned values. On 
the first call, T is set equal to TZRO, TEND is set equal to TZRO + DTEND, and 
the Y array is assigned the initial values for the dependent variables. The 
initial value of the step size, DELH, also must be chosen. This value should 
be considerably smaller than the average value expected for the problem since 
DGEAR must begin with a not so highly efficient first-order method. For the 
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first step, as for every step, DGEAR tests for the possibility that the step 
size was too large to pass the error test (based on TOL) , and so adjusts the 
step size down automatically. This downward adjustment, if any, is noted by 
IER having the values of 66 or 67. The relative error bound, TOL, should be 
assigned a value, TOL is an estimate of the local truncation error, as dis- 
tinguished from the global truncation error. The local truncation error is 
the error resulting from taking a single step, starting with data which is 
regarded as exact. TOL should be at least an order of magnitude larger than 
the unit roundoff but generally not larger than 0.001. 

DGEAR solves the differential equations from T = TZRO to T = TZRO 4- DTEND 
with the given initial conditions. The value of IER, which is returned to the 
MAIN program from DGEAR, is tested and the program stops if IER is between 132 
and 136. The output error parameter, IER, informs the user of errors encount- 
ered by the GEAR package, if any. If the program failed to pass the error test 
after reducing DELH by a factor of 1.0E10 from its initial value, then IER is 
set equal to 132. If the program failed to achieve corrector convergence after 
reducing DELH by a factor of 1.0E10 from its initial value, then IER is set 
equal to 133. Incorrect input to the GEAR package is flagged by IER = 135 or 
136. If IER is less than 132 or greater than 136, the desired output is 
printed, TEND is reset to TEND = TEND + DTEND, and DGEAR is called again. This 
process is continued until there is seal failure, TEND = TSTOP, or IER is be- 
tween 132 and 136. To solve the problem, the GEAR package used the implicit 
Adams Method and used the following values for its parameters; N6 = 6, 

DELH = 0.001, TOL = 0.0001, METH = 1, MITER = 0, and INDEX = 1. 

To find the moments M^ and M 2 and the force F z , the double integrals 
of Eqs. (8), (9), and (10) have to be evaluated. This was done by the applica- 
tion of the Guassian quadrature integration formula, one of a family of powerful, 
accurate, numerical integration techniques £ 3 ^. To evaluate these integrals, 
the Fortran program SQUAD1 was used. This program was written at the Lewis 
Research Center by Charles Goldstein, and it evaluates integrals by repeated 
application of Gauss 1 formula 01 

The call to obtain the inner integrals, which give integration with respect 
to R, is CALL SQUADl (MODE, NR, RI, DDl, XR, YR1, ANSI), where the parameters 
are defined as follows. 

MODE - =1, then obtain the arguments XR(NR) 

=2, then use the evaluations of the integrand YRl[XR(NRj] to 
obtain ANSI 

NR - The number of points of integration; N = 3(1)16 

RI - The lower limit of integration 

DDl - The upper limit of integration 

SR - An array containing the variable of integration 

YR1 - An array containing values of the integrand 
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ANSI - The result of the integration; i.e., the value of the integral 

The call to obtain the outer integrals, which gives integration with re- 
spect to 0, is CALL SQUAD1 (MODE, NTH, THl, TH2, XTH, YTH1, ANSI), where the 

parameters are defined as follows. 

MODE - =1, then obtain the arguments XTH (NTH) 

=2, then use the evaluations of the integrand YTHl[XTH(NTH)]] to 
obtain ANSI 

NTH - The number of points of integration per quadrant; NTH = 3(1)16 

THl - The lower limit of integration 

TH2 - The upper limit of integration 

XTH - An array containing the variable of integration 

YTHl - An array containing values of the integrand 

ANSI - The result of the integration; i.e., the value of the integral 

The program used 5 points of integration for the inner integral (NR = 5), and 
8 points of integration for each quadrant of the outer integral (NTH = 8). 

Larger values were tried for NR and NTH (e.g., NR = 11 and NTH = 16) and the 
values of the resulting integrals were compared. It was found that using 
NR = 5 and NTH = 8 gave sufficient accuracy for the integrals. The variables 
NR and NTH can take on values between 3 and 16. To compute the outer integrals, 
whose limits of integration are 0° and 360°, SQUADl was used four times over the 
quadrants 0° to 90°, 90° to 180°, 180° to 270°, and 270° to 360°. This was done 
because NTH can have a maximum value of only 16, and applying SQUADl over four 
intervals allowed more points to be used over the total interval of integration 
0° to 360°. It was also done so that the program would avoid computing the 
integrand at the points 0° and 180°. The routine SQUADl does not compute the 
integrand at the end points of the interval over which SQUADl is being used, so 
that computation at these points was avoided. Singularities can only occur at 
0° and 180°, but because of the complexity of the integrand, it could not easily 
be determined whether singularities do occur. As a precaution, computation of 
the integrand was thus avoided at 6=0° and 180°. 


FORTRAN PROGRAM 
General Description 

The foregoing analysis has resulted in a computer program for the analysis 
of a noncontacting mechanical face seal. The MAIN executive program handles 
most of the program flow and logic. It assigns values to program parameters, 
reads and writes program input, initializes the GEAR package and repeatedly calls 
the GEAR package to move ahead in time by DTEND at each call. After each call, 
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it writes the output and tests for seal failure or a failure of the GEAR package. 
The program stops when a failure is detected or when TSTOP is reached. The sub- 
routine FQUAD does most of the computing for the functions which are used in the 
DGEAR package. 

A brief explanation of the function of each program module follows. 

DGEAR the general GEAR package for solving a system of first order 
ordinary differential equations 

ERMIN a subroutine called by ERTST to determine when Hw varies by 

less than HSTAR during several revolutions. If this condition 
is met, then ERMIN sends a flag through ERTST to MAIN 

ERTST a subroutine called by the MAIN program to compute H M j n and 
determine if there is seal failure or GEAR package failure 

FCN a subroutine which computes the Y array 

FCNJ a dummy subroutine passed to the GEAR package 

FQUAD a subroutine used to compute the integrand of the double integrals 

appearing in Eqs. ( 8 ) through (10) 

MAIN the Main executive program 

SQUAD1 a subroutine called by FQUAD to perform a single Gaussian 
integration 

WRT a subroutine called by the MAIN program which writes the output 

Y 9 ip 9 Z, y, 4 * 9 > a y > HFiXN ^ T, ^x> a x> a y* fty> Z, DELH, 

INDEX, IER, and FC after each call to the GEAR package 

Using the Program 

The program always starts by reading program input through the two name- 
lists, INPUT1 and INPUT2. The input in these namelists are described below (all 
symbols are repeated in appendix C) . 

Namelist INPUT1 


DELH a parameter in the GEAR package containing the first step size to 
be tried 

DTEND for each call to the GEAR package, solve the differential equa- 
tions from TEND to TEND + DTEND 

HMSTOP the program is stopped when HMIN becomes less than HMSTOP 

HSTAB the program is stopped if HMIN vairies by less than HSTAB for 

several revolutions 
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INDEX a parameter in the GEAR package indicating the type of call being 
made to the package 

METH a parameter in the GEAR package call indicating the method of 

solution to be used 

MITER a parameter in the GEAR package call indicating the iteration 
method to be used 

NPLT if punch is true, then write every NPLT-TH array point into data- 

set 15 to be saved to make a movie 

NR a parameter in the SQUAD1 call indicating the number of points of 

integration used to obtain the inner integrals 

NTH a parameter in the SQUAD1 call indicating the number of points of 

integration used to obtain the outer integrals 

PUNCH a logical variable. If PUNCH = TRUE, then write and save output 
in dataset 15 

READ1 a logical variable. If READ1 = TRUE, then read in and process 
another case 

TOL a parameter in the GEAR package call giving the relative error 

bound 

TSTOP the differential equations are solved to TEND = TSTOP 

TZRO the initial value of time, the independent variable 

YIN an input array containing the initial values for ct x , a y , Z, cf x , 

Gty , and Z 

Namelist INPUT 2 

BETA normalized coning angle, 3 

Q 

FKROC the dimensionless parameter, K which is the spring constant 

L o 

multiplied by the equilibrium center-line clearance divided by r Q 

FM dimensionless mass, m 

FME dimensionless external moment, M e 

FI dimensionless pressure at the inner radius of the ring, Pj_ 

FO dimensionless pressure at the outer radius of the ring, P Q 

RG dimensionless ring radius of gyration, Rg 



RI 

RSP 


11 

dimensionless inner radius, R^ 

dimensionless radius of the location of the springs, R S p 


If the logical variable READ1 is TRUE, then another seal case is read into 
the MAIN program through the Namelists INPUT1 and INPUT2 and this case is pro- 
cessed. If READl is read in as FALSE, then the program stops. 
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APPENDIX A 
NOMENCLATURE 

(1 - R) /H m H 2 (l - R.) 
seal center-line clearance 
equilibrium center-line clearance 
closing force 

dimensionless closing force, F*/Sr^ 
axial force 

dimensionless axial force, F*/Sr2 
dimensionless film thickness, h/C 
dimensionless minimum film thickness 
film thickness 

ring mass moment of inertia about a diameter 
dimensionless moment of inertia, I* 0 )^C o /Sr^ 
spring constant 

dimensionless spring constant, K*/Sr Q 
moment 

external moment 

dimensionless external moment 

dimensionless moment, M*/Sr~^ 

o 

ring mass 

dimensionless mass, m*m^C /Sr^ 

o' o 

number of seal revolutions, t / 2 tt 
dimensionless pressure, p/S 
pressure 

dimensionless radius, r/r Q 
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r radial coordinate 

r g ring radius of gyration 

Rg dimensionless ring radius of gyration 

5 seal parameter, 6yu(r 0 /C 0 ) (1 - 

t * t ime 

t dimensionless time, oat 

Z* axial displacement 

Z dimensionless displacement, Z X /C Q 

a* tilt angle 

cl normalized tilt, a r 0 /C 0 

3 coning angle 

8 normalized coning, 8 r 0 /C 0 

y nutation 

Y normalized nutation, y r 0 /C 0 

6 coning parameter, g*r Q /C 

e tilt parameter, y*r 0 /C 

8 angular coordinate 

y viscosity 

^ precession 

oj shaft angular velocity 

Subscripts : 

1,2,3 axes 1, 2, or 3, respectively 
d hydrodynamic 

i inner radius 

m mid radius 


o 


outer radius 
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s 

sp 

x,y,z 


hydrostatic 

springs 

axes, x, y, or z, respectively 
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APPENDIX B 
FORTRAN SYMBOLS 

ANSI a parameter in the SQUAD1 call containing the value of the integral 
BETA normalized coning angle, B 


CAX 

normalized tilt , a , 

X 

wrt 

. the 

x axis 

CAY 

normalized tilt, a^, 

wrt 

. the 

y axis 

CAXD 

the first derivative 

of 

a 

X 


CAYD 

the first derivative 

of 

a y 



CZ the axial dimensionless displacement of the ring, Z 

CZD the first derivative of Z 

DD1 the upper limit of integration in the SQUAD1 call 

DEL the coning parameter, 6 

DELH a parameter in the GEAR package containing the step size to be used on 
the first call to DGEAR and the step size actually used on subse- 
quent calls 

DTEND for each call to the GEAR package, solve the differential equation 
from TEND to TEND + DTEND 

EPS the tilt parameter, e 

FC the dimensionless closing force, Fc 

FCN a subroutine used by the GEAR package for evaluating the function Y 

FCNJ a DUMY subroutine used by the GEAR package 

FI the dimensionless moment of inertia, I 

C 0 

FKROC the parameter K — , which is the dimensionless spring constant multi- 

r o 

plied by the equilibrium center-line clearance divided by r Q 
FM dimensionless mass, m 

FM1 dimensionless moment, M^, wrt. axis 1 

FM2 dimensionless moment, M 2 , wrt. axis 2 

FME dimensionless external moment, M e 
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FMX 

FMY 

FZ 

GAM 

GAMD 

H 

HI 

HO 

HM 

HMIN 

HMIN1 

HMIN2 

HMSTOP 

HSTAB 

IER 

INDEX 

IWK 

IWRT 

LOGER 

METH 

MITER 


dimensionless moment, wrt. axis x 
dimensionless moment, My, vrrt. axis y 
dimensionless axial force, F z 
normalized nutation, y 

the first derivative of the normalized nutation, y 
dimensionless film thickness, H 

dimensionless film thickness at the inner radius, H^ 
dimensionless film thickness at the outer radius, H Q 
the average of HI and HO 

dimensionless minimum film thickness, H^in 

dimensionless minimum film thickness at the inner radius, H^^]_ 

dimensionless minimum film thickness at the outer radius, H^-j-^ 

the program is stopped when HMIN becomes less than HMSTOP 

the program is stopped if HMIN varies by less than HSTAB for several 
revolutions 

an output from the GEAR package indicating any errors encountered 
while using the GEAR package 

an output parameter from the GEAR package call indicating whether the 
integration was successful, and an input to the GEAR package indi- 
cating the type of call being made to the package 

a parameter in the GEAR package call which is used as an integer work 
vector 

a parameter indicating whether debug output is wanted 

a logical variable returned from subroutine ERTST. LOGER = FALSE means 
that the program should be stopped 

a parameter in the GEAR package call indicating whether the Adams 

method or the backward differentiation formulae method is to be used 

an input in the GEAR package call indicating whether the chord method 
or functional iteration is to be used and whether a finite difference 
or analytic Jacobian is to be used with the chord method 
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N6 a parameter in the GEAR package call indicating the number of first 

order differential equations to be solved 

NPLT a program input indicating how much output to save for movie 
generation 

NR a parameter in the SQUAD1 call indicating the number of points of 

integration used to obtain the inner integrals 

NTH a parameter in the SQUAD1 call indicating the number of points of 

integration used to obtain the outer integrals 

P dimensionless pressure 

PD the dimensionless hydrodynamic and squeeze pressure components, P^ 

PI the dimensionless pressure at the inner radius of the ring, P^ 

PO the dimensionless pressure at the outer radius of the ring, P Q 

PS the dimensionless hydrostatic pressure component, P s 

PSI the precession, ip 

PSID the first derivative of the precession, ip 

PUNCH a logical variable. If punch = true, then write output into dataset 15 

to be saved to make movies 

RAX, RAY, 

RTANG, 

RZ 

READ1 

RG 
RI 
RM 
RSP 
T 

TEND 


arrays in which to save output for movie generation 

a logical variable. If READl = TRUE, then read in and process another 
case 

the dimensionless ring radius of gyration 

the dimensionless inner radius, R^ 

the dimensionless mid radius, R m 

the dimensionless radius of the location of the springs, R 

sp 

dimensionless time, the independent variable 

a parameter _to the GEAR package giving the value of T at which a solu- 
tion for Y is desired 

a parameter in the SQUAD1 call giving the lower limit of integration of 
the outer integral 


TH1 
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TH2 a parameter in the SQUAD1 call giving the upper limit of integration 

of the outer integral 

TOL a parameter in the GEAR package call giving the relative error bound 

TOUCH a logical variable. TOUCH = T means that the seal ring touches the 

shaft 

TSTOP the differential equations are solved to TEND = TSTOP 

TZRO the initial value of time, the independent variable 

WK a parameter in the GEAR package call containing a real work vector 

XR a parameter in the SQUAD1 call containing the variable of integration 

of the inner integrals 

XTH a parameter in the SQUAD1 call containing the variable of integration 

of the outer integrals 

Y a parameter in the GEAR package call containing the dependent variable, 

Y 

YIN an input array containing the initial values for a x , ay, Z, a # x , a y , 

and Z 

YPRIME a parameter in the FCN call containing the first derivative of the Y 
array 

YR1 a parameter in the SQUADl call containing the values of the integrand 

of the inner integrals 

YTHl a parameter in the SQUADl call containing the values of the integrand 

of the outer integrals 
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* 


APPENDIX C 

FORTRAN PROGRAMS 


1 JC 

2 :c 
3:c 
4tc 
s:c 
6 : c 
7:c 
8:c 
9:c 
io:c 
n :c 
i 2 :c 
i3:c 
1 4 t C 

is : c 
i6:c 
17 : c 
i8:c 

19JC 

2 o:c 
2 1 : c 
22 :c 

23JC 

24JC 

25JC 

26 : c 

27tC 

28 :c 

29 :c 
3o;c 

31 :c 

32 :c 
33JC 
34 :c 
3s:c 

36 1 C 

37 t C 
38:C 
39 J C 
4o:c 
4i:c 
42JC 
43:c 
44:c 
45JC 
46:c 
47JC 
48tC 
49 J C 

so:c 

51 :c 

52 ! C 
535C 
54 5 C 
55:c 
56:c 

57 J C 

58 :c 
59:c 
6o:c 
6i :c 


MAIN PROGRAM 

THIS PROGRAM DESCRIBES THE MOTION OF A SEAL RING 

BY SOLVING SIX FIRST ORDER DIFFERENTIAL EQUAS ♦ (WHICH 
WERE DERIVED FROM THREE 2ND ORDER ORDINARY DIFFERENTIAL 
EQUATIONS) 


THE SYSTEM OF DIFFERENTIAL EQUAS, 
YPRIME(l) = AX -DERI 
YPRI ME ( 2 ) = AY-DER1 
YPRIME<3) = Z-DER1 


YPRIMEC4) 
YPRIME(5) 
YPRIME ( 6 ) 


= AX-DER2 
= AY-DER2 
= Z-DER2 


TO BE SOLVED IS AS FOLLOWS 


WHERE THE Y ARRAY IS DEFINED AS FOLLOWS 

Y < 1 ) =AX 

Y ( 2 ) = A Y 
Y ( 3 ) =Z 

Y ( 4 ) = DER-AX 
Y ( 5 ) = DER-AY 
Y ( 6 ) = DER-Z 


DTEND INCREMENT TEND BY DTEND 


HSTAB - STOP THE PROGRAM IF HMIN VARIES LESS THAN HSTAB 
DURING FIVE REVOLUTIONS 

HMSTOP - STOP THE PROGRAM IF HMIN BECOMES LESS 
THAN HMSTOP 


I WRT = 0 DO NOT WRITE DEBUG OUTPUT 

= 1 WRITE COMPLETE DEBUG OUTPUT 


NPLT - PUNCH EVERY NPLT FT, FOR THE MOIVE 
PLOTTING PROGRAM 

NR = NO* OF INTEGRATION INTERVALS IN THE 
INTERVAL 1 TO RI 


NTH = NO. OF INTEGRATION INTERVALS 

IN EACH QUADRANT OF THE INTERVAL 0 TO 2*PIE 

PUNCH - PUNCH= ♦ T * * THEN PUNCH OUTPUT FOR THE MOVIE 
PLOTTING PROGRAM 

PUNCH= ♦ F ♦ f THEN DO NOT PUNCH THIS OUTPUT 
FOR THE MOVIE PLOTTING PROGRAM 


RAX>RAY>RZ»RTANG - INFORMATION PUNCHED FROM THIS PROG* 
TO BE READ INTO l USED IN THE MOVIE PLOTTING 
PROGRAM 

TEND FOR EACH CALL TO DGEARf FIND SOLUTION 
AT TEND 

TOUCH= * T « MEANS THAT THE SEAL RING TOUCHES THE SHAFT. 


TSTOP SOLVE DIFF. EQUAS* TO T=TSTOP 
TZRO INITIAL T 
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62;c 
63 :c 

64JC 

65;c 

66: 

67:c 

68 : 

695C 

70: 

71 : 

72: 

73 JC 

74: 

75 : c 

76: 

77: 

78:c 

79: 

80:c 

site 

82JC 
83: 
84: 
85: 
865 
87:c 
88: 
89 i 
9o:c 
91 : 
92JC 
93: 
94 :c 
95JC 
96: 
97:c 
98:c 
99: 

100 : 
101 : 
102: 
103: 

104 1 

105: 

106: 

107: 

108: 

io9:c 

no: 

111 : 

112: 

113: 

114: 

115: 

116: 

117! 

ns: 

119: 

i2o:c 

121 :c 

122 : 

123: 

124: 

125: 

126: 

127: 

128: 


EXTERNAL FCN * FCN J 

LOGICAL LOGERf TOUCH f READ! » PUNCH 

DIMENSION Y I N ( 6 ) 

DIMENSION Y ( 6 ) f UK ( 1 44 ) f IUK ( 6 ) 

DIMENSION RAX ( 6200 > f RAY (6200) f RZ ( 6200 ) r RTANG ( 6200 ) 

COMMON/COM 1/ PIE f F*IE2 ? N4PIE f NR » NTH » IURT 

C0MM0N/C0M2/ FKROC r RO f BETA f R I t CO » 

1 RMfPOf PI fFI f FMEfFMfRSP 

C0MM0N/C0M3/ HMIN t HMIN1 f HMIN2 

COMMON/GEAR/ DUMMY (48) fSDUMMY(4) * I DUMMY (38) 

NAMELIST /IN PUT1/ READ 1 * METH * M I TER * TOL » I NDEX f DELH f 

1 TZROf YINf TSTOPr 

2 DTENDf NR 1 NTH* HMSTOP f HSTAB f 

3 NPLT t PUNCH 

NAMELIST /INPUT 2/ RI f FKROC t BETA » PI f PO r 
1 FMEfFMfRGf RSP 

DATA RE ADI / . TRUE ♦ / 

DATA TZRQfTSTOPt DTENDf YIN /O ♦ f 1 3 ♦ f , 1 f 0 . f . 2 f 4*0 , / 


DATA HMSTOP/. IE-20/ 


NR = 5 
N 4 F‘ I E = 4 
NTH = 8 
N6 = 6 
I URT = 0 

PUNCH® « FALSE ♦ 

NPLT = 1 

PIE=3. 14159265359 
PIE2=2«#PIE 
HSTAB® . 00 1 

GO TO 14 

20 I F ( ♦ NOT * PUNCH ) GO TO 14 

WRITE ( 15 f 122) RI t FKROC t BETA f PI fPOf FMEfFM 

WRITE ( 15 f 124 ) RGfRSPfTZROfTOLfDELHfNRfNTHfMETHf MITER 

WRITE(15f 122) YIN 

WRITE(15f 112) JTfNPLTfTOUCH 

WRITE ( 15 f 114) ( R A X ( J ) f J= 1 f JT f NPLT ) fRAX( JT) 

WRITE < 15 f 114) ( RAY ( J ) f J= 1 f JT f NPLT ) fRAY( JT) 

WRITE ( 15 f 114) (RZ( J) f J= 1 f JT » NPLT ) fRZ( JT) 

WRITE ( 15 » 114) (RTANG ( J ) f J= 1 f JT f NPLT ) f RTANG ( JT ) 


14 METH= 1 
M I TER=0 
INDEX® 1 
TOL®. 0001 
DELH® .01 

I F < ♦ NOT • READ1 > GO TO 99 
REAI»( 5 f INPUT 1 ) 
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129: W R I T E ( 6 * 101 ) 

130 J URITEUfINPUTI) 

131 J READ ( 5 » INPUT2 ) 

132: UIRITE(6fI NPUT2 ) 

1 33 : c 

134JC 

135J FI=*5*FM*RG*RG 

136: RM = ( 1 • +RI ) /2 ♦ 

1 37 : c 

138:c 

139:c INITIALIZE THE INTEGRATION 

i4o:c 

1 4 1 J 30 LOGE R= * TRUE ♦ 


142: 

143: 

1 44 : c 

i4s: 
146: 
147 JC 
148: 
149: 
150 : 
151 : c 
i52:c 
153IC 
154: 
155: 
156: 
i57:c 
158: 
i59;c 
160: 

1 6 i : c 

162 J 

163 J 
164 : 
165: 
166: 

1 67 : c 

1 68 : c 

169IC 
170 : 
171 : 
172: 

1 73 ; c 

174: 
175: 
176: 
177: 
178: 
i79:c 
180 1 
i8i:c 
182: 
i83:c 
184 J 
185:c 
186:c 

187:c 
IBs: 
189:c 
190: 
191 : 

192JC 

193:c 

194:c 

195: 


TOUCH=. FALSE* 
JT= 1 

T = T Z R 0 

TEN D=TZRO+D TEND 

DO 32 J32=l f 6 
Y (J32)=YIN(J32) 
32 CONTINUE 


WRITE OUTPUT HEADING 
WRITE(6f 104) 

URITEC7, 106) 

UR I TE ( 7 f 102 ) 

CALL ERTST < LOGER f TOUCH f Y f IERfT f DELH f DTEND f HMSTOP f HSTAB > 

CALL WRT(N6f JTt Yf Tt DELHi TEND f INDEX f I ER) 

IF ( ♦NOT. PUNCH) GO TO 40 
R AX ( JT > = Y ( 1 ) 

RAY ( JT )=Y ( 2) 

R Z ( J T ) = Y ( 3 ) 

RTANG(JT)=TZRO 

USE SUB DVOGER TO MOVE ONE STEP IN TIME 
40 JT= JT+ 1 

CALL DGEAR(N6fFCNfFCNJfTfDELHfYfTENDfT0Lf 
1 METHf MITER f INDEX f IUKfUKf IER) 

I F ( .NOT. PUNCH) GO TO 54 
RAX ( JT ) =Y < 1 ) 

RA Y ( JT ) = Y ( 2 ) 

RZ ( JT ) = Y ( 3 ) 

RTANG( JT)=TEND 

54 CAtL ERTST(LOGERfTOUCHf Yf IERf T fDELH f DTENDf HMSTOP f HSTAB > 
I F ( .NOT. LOGER) GO TO 20 

CALL WRT(N6f JTf YfTf DELHf TEND t INDEX f IER) 


TEN D=TEND+D TEND 

IF(TEND.GE.TSTOP) GO TO 20 
GO TO 40 


101 FORMAT ( 1H1 ) 
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196: 104 FORMAT ( 1H1 f 6Xf 'C Z ' f 9X f ' GAM ' f 8X f ' F’SI ' f 8X f 'CZD' fBXf 

197 5 1 'GAMD' f 7 X f ' P S I D ' »7X? 'CAX ' f8Xf 'CAY 7 f 8X f ' HMI N ' f 7X f 

198J 2 ' TWRT ' ) 

1 99 : c 

200 : 102 FORMATCIHIf 'INDEX 'f2Xf 'IER'fSXf 'AXD'f8Xf 'AYD'f9Xf 

201 : 1 'ZD' f9Xf 'T' f 10X f 'H' ? 9X f ' AXDD ' f 7Xf ' AYDD' f7Xf 

202: 2 ' ZDD ' f 9X f ' FC ' ) 

203:C 

204: 106 F0RMAT(1H0f/////) 

205: 112 FORMAT ( 2 1 4 f L4 ) 

206: 114 FORMAT (8E10.4) 

207: 122 FORMAT ( 7E1 1 * 4 ) 

208: 124 FORMAT (5E11*4f4I5) 

209:C 

210: 99 STOP 

21 1: END 
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14 
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19 

20 
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23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 
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56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 


SUBROUTINE ERTST (LOGER > TOUCH f Y , IER i T i DELH > DTEND > HMSTGP > HSTAB ) 


THIS SUBROUTINE TESTS THE VARIABLE IER 
FROM THE DGEAR CALL AND PRINTS 
THE CORRECT ERROR MESSAGE 

THE SUBROUTINE ALSO COMPUTES HMIN 


DIMENSION Y < 6 ) 

LOGICAL LOGER » TOUCH 

COMMON/COM 1 / PIE » PIE2 * N4PIE? NR *NTH » IWRT 

C0MM0N/C0M2/ FKROC»RO*BETA»RI>CO» 

1 RMrPOfPIfFlrFMEfFMfRSP 

C0MM0N/C0M3/ HMINtHMINl , HMIN2 

TO COMPUTE AND TEST HMIN 


CAX = Y ( 1 ) 

C A Y=Y ( 2 ) 

SQXY=CAX*CAX+CAY*CAY 
G AM= SORT ( SQXY ) 

C AXX= ABSCCAX) 

C A Y Y= ABS(CAY) 

PSI= ATAN2(CAYYtCAXX) 

IF ( CAX ♦ LT ♦0*AND.CAY.GT *0) PSI=-PSI 
I F < C AX . GE • 0 ♦ AND « CAY • LT « 0 ) PSI = -PSI 

CPSI=CAX/GAM 

CZ1=1 .+Y<3> 

EPS=GAM/CZ1 
DEL=BETA/CZ 1 

HMIN1=1 .-EPS*RI 
HMIN2=1 *-EPS+DEL#(l*-RI ) 

HMIN=HMIN1 

IF (HMIN2.LT ♦ HMIN1 > HMIN=HMIN2 
IF ( HMIN ♦ GT *HMSTOP) GO TO 20 
L0GER= ♦ FALSE ♦ 

TOUCH= ♦ TRUE ♦ 

WRITE(6f 102) HMIN* HM INI ? HMIN2 
GO TO 99 

20 HI =HMI N 

CALL ERMI N< DTEND, HI f LOGER » HSTAB) 

22 IF ( IER.NE.33) GO TO 24 
LOGER= . FALSE ♦ 

WR I TE < 6 , 104) IER 
GO TO 99 

24 IF ( I ER ♦ LT « 132«0R.IER*GT » 136) GO TO 99 
LOGER=. FALSE. 

UR I TE ( 6 , 104) IER 


101 FORMAT ( 1H1 ) 

102 FORMAT ( ' FATAL ERROR ' , 5X , / HMIN= / E12.5»5Xf 
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68J 1 'HMINl='E12.5f5Xf 'HMIN2='E12.5) 

6?i 104 FORMAT (' FATAL ERROR '» 5X r ' IER= ' 15) 

70:c 

71 J 99 RETURN 
72: END 



25 


1 J SUBROUTINE ERMIN ( DTEND r Hi f L06ER f HSTAB ) 

2JC 

3:c 

4JC THIS SUBROUTINE TESTS HNIN ♦ IF HhIN DOES 

5 J C IN FIVE REVOLUTIONSf THEN SET 

6JC 

7:C IN CALLING PROG, ERTST t HI -HNIN 

8:c 

9JC 

io:c 

1 1 : DATA IDOfHS /Of 1 « E20/ 

1 2 J C 

13J LOGICAL LOGER 

1 4 : c 

15: IF ( IDO « EG « 1 ) GO TO 40 

i6:c 

17J IF( ABS(Hl-HS) *GT. HSTAB) GO TO 50 

18 :c 

19: HS=H 1 

20: TS=0 ♦ 

21 : ID0=1 

22: GO TO 99 

23:c 

24 : so i d o = o 

25: HS=H 1 

26 J GO TO 99 

27 : C 

28: 40 IF ( ABS ( HS-H1 ) ♦ LT .HSTAB ) GO TO 55 

29: I D 0 = 0 

30: HS=H 1 

31 : GO TO 99 

32:c 

33: 55 TS=TS+ DTEND 

34: IF ( TS ♦ GT ♦ 3 1 ♦ 4 16 ) LOGER= * FALSE ♦ 

35 :c 

36: 99 RETURN 

37: END 


NOT VARY BY HSTAB 
LOGER TO ♦ FALSE « 
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I! SUBROUTINE FCN< N t T iYi YF'RIME ) 

2 5 C 

3 5 C THIS SUB IS A PARAMETER IN THE CALL 

4 SC TO HGEAR 

5SC 

65 C THIS SUB COMPUTES YF'RIME 

7 * C 

85 DIMENSION Y(N) » YPRIME(N) 

9 • C 

105 C0MM0N/C0M1/ PIE , PIE2 , N4PI E . NR » NTH , IWRT 

1 1 : c 

125 C0MM0N/C0M2/ FKROC i RO » BET A » RI » CO » 

135 1 RH»PO»PI,FI,FMEfFM,RSP 

1 4 5 C 

155 C0MM0N/C0M4/ CAX * CA Y , CZ » CAXD> CA YD t CZD 

165 COMMON/COH5/ FMX » FMY » FZ t GAM » PSI » GAMD , PSID » FC f 

175 1 FM1 »FM1C»FZ1»FM2» EPS>DEL 

1 8 5 C 

19 5 C 

20 5 C COMPUTE THE RIGHT HAND SIDE OF THE 6 FIRST ORDER 

21 5 C DIFFERENTIAL EQUATIONS 

22 5 C 

235 20 CAX= Y ( 1 ) 

245 CAY=Y ( 2 ) 

255 CZ = Y ( 3 ) 

265 CAXD = Y ( 4 ) 

275 CAYD=Y ( 5 ) 

28 5 CZD = Y < 6 ) 

29 5 C 

305 CALL FQUAD 

31 5 C 

325 YPRIME ( 1 ) =Y< 4 ) 

33 5 YF'RIME ( 2 ) =Y < 5 ) 

34 5 YF’RI ME ( 3 ) = Y ( 6 ) 

355 YF‘RIME(4) = (FMX+FME)/FI 

365 YPRIME(5)=FMY/FI 

375 YPRIME(6)=FZ/FM 

38 5 C 

39 5 C 

405 10 RETURN 

415 END 
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i: SUBROUTINE FCN J ( N > T f Y , PD ) 

2tC TO COMPUTE THE N*N JACOBIAN MATRIX OF PARTIAL 

3tC DERIVATES FOR THE DGEAR IMSL ROUTINE 

4JC THIS IS A DUMY ROUTINE FOR DGEAR 

5JC 

6: DIMENSION Y(N)*PD(NiN> 

7JC DOUBLE PRECISION T»Y»PD 

8:c 

9: DUMY=DUMY 

10: RETURN 

11 : END 



28 


1 J 

2 ♦ C 

3 ♦ C 

4 : c 
5:c 
6 t C 
7:c 

8:c 

9:c 

io:c 

u:c 

i 2 :c 

13JC 
1 4 t C 

is:c 

i6:c 

i7:c 

18 * C 

i9:c 

20 : 

21 : 

22JC 
23: 
24JC 
25: 
26: 
27:c 
28: 
29 : c 

30 : 

31 : 
32:c 
33:c 
34: 


SUBROUTINE FQUAD 


THIS SUBROUTINE PERFORMES DOUBLE INTEGRATION 
USING THE GOLDSTEIN QUADRATURE ROUTINES 
TO FIND Ml*M2rAND FZ 

INTEGRATE THE INNER INTEGRAL WRT R 
INTEGRATE THE OUTER INTEGRAL WRT THETA 


NR = THE NO* OF PTS. IN THE R INTEGRAL 
NTH = THE NO. OF PTS IN EACH QUADRANT OF THE THETA 
INTEGRAL. WE INTEGRATE FROM 0 TO 2*PIEr 
& USE NTH PTS. IN EACH QUADRANT 

I WRT = 0 DO NOT WRITE DEBUG OUTPUT 
=1 WRITE DEBUG OUTPUT 


DIMENSION XTH( 16) > XR < 1 6 ) » YTH1 (16) ? YTH2 (16) f 
1 YTH3( 16) r YR1 ( 16) v YR2< 16 > v YR3U6) 

COMMON/COM 1/ PIE > PIE2 f N4PIE f NR f NTH # I WRT 

CDMM0N/C0M2/ FKROC , RO , BET A > R I * CO f 
1 RMfPOfPIfFIrFMErFhfRSP 

C0MM0N/C0M4/ CAX f CA Y » CZ > C AXD f C AYD f CZD 

COMMON/COM5/ FMX f FMY f FZ * GAM r PS I f GAMD f PS I D f FC f 
1 Fhl f FM1 C f FZ 1 f FM2 f EPS f DEL 


FN4FIE=N4PIE 


35J PIEN4=2.*PIE/FN4PIE 

36:C 

37*. D D 1 = 1 ♦ 

38: Fh 1 =0 « 

39: FM2=0 ♦ 

40: FZ 1 =0 ♦ 

41 :c 

421 SQXY=CAX*CAX + CAY*CAY 

43J GAM= SQRT(SQXY) 

44JC 

451 CAXX= ABS(CAX) 


46J CA YY= ABS(CAY) 

47: PSI= ATAN2(CAYY»CAXX) 

48: IF(CAX.LT.O.AND.CAY.GT.O) PSI=-PSI 

49: IF (CAX.GE.O.AND.CAY.LT *0) PSI=-PSI 

so:c 

51 : CPSI=CAX/GAM 

52: SPS I =C AY/GAM 

53:c 

54: GAMD=CAXD*CPSI + CAYD*SP5I 


55: PSID=(CAYD*CPSI - CAXD*SPS I ) /GAM 

56: IF(IURT.GE.l) WRITE ( 1 1 f 1 16 ) C AX f CAY f CZ f CAXD f CAYD f CZD f GAM f GAMD 

57:c 

58: RI Ml = 1 ♦ -RI 

59: Z1=1.+CZ 

60: EPS=GAM/Z1 

61 : DEL=BET A/Zl 

62: EPSRI=EPS*RI 

63: DELRI=DEL*RIM1 

64: Z3R= 1 ./Z1/Z1/Z1/RIM1 

65J PD 1 = ( .5-PSID)*GAM*RM 

66 4 . PD2=GAMD*RM 

67J PS 1 = • 5* ( PO-PI ) 
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CALL SQUAD 1 (IfNRfRIfDDIfXRfYRIfANSI) 

TH1 =0 ♦ 

TH2=0. 

I F ( I WRT « EQ « 1 ) WRITE ( 1 1 f 102 ) NR f RI f ANSI f PSI f CPSI f SF’SI f 

1 GAM, GAhDrPSIDfZl rEPSf DELfEPSRIfDELRIfZ3R»PDl»F'D2f 

2 PSI fXRfYRl 

DO 40 J4 = 1fN4F‘IE 

TH1=TH2 

TH2=TH2+PIEN4 

CALL SQUAD1 (1 fNTHfTH1fTH2fXTHfYTH1fANS1) 

IF(IURT.EQ.l) WRITE ( 1 1 f 1 04 ) J4 f NTH , TH1 f TH2 f ANSI f XTHf YTH 1 
DO 30 J3= 1 f NTH 
THA=XTH ( J3 ) 

STHA = SIN(THA) 

CTHA= COS(THA) 

DO 20 J2=l f NR 
R=XR ( J2 ) 

R2=R*R 

ECTHA=EPS*CTHA 

H=1 .+R*ECTHA+DEL*(R-RI ) 

H I = 1 .+EF'SRI*CTHA 
H0 = 1 . +ECTHA+DELRI 
HM = ♦ 5* < H I +H0 ) 

A = < 1 ♦-R)/<HM*H*H*RIM1 ) 

PD=Z3R*(PD1*STHA“CZD-PD2*CTHA>*A*(R-RI ) 

PS=P0-PS1*HI*HI*(H0+H)*A 

F‘=PD + PS 

I F ( P * LE ♦ 0 ) F’ = 0 * 

YR1C J2)=P*R2*CTHA 
YR2< J2)=P*R2*STHA 
YR3< J2)=P*R 

IF ( I WRT • EQ • 1 ) WRITE ( 1 1 f 106) R f R2 f ECTHA f H f HI f 
1 HGfHMf Af PDfPSfP 
20 CONTINUE 

CALL SQUAD1 <2fNRfRIfDD1f XRf YR1 fANSI ) 

CALL SQUAD1(2f NRfRIfDD1fXRfYR2fANS2) 

CALL SGUAD1 <2fNRfRI fDDI fXRfYR3f ANS3) 

YTH1 ( J3 ) = ANSI 
YTH2 ( J3 ) =ANS2 
YTH3 ( J3 ) =ANS3 


IF ( I WRT » EQ « 1 ) WRITEdlf 108) NTH f THA * STHA f CTHA f ANSI , 
1 ANS2f ANS3fYR1 fYR2fYR3 

30 CONTINUE 

CALL SQUAD1 ( 2 f NTH f TH1 f TH2 » XTH f YTH1 f ANSI) 

CALL SQUAD1 (2fNTHfTH1 fTH2fXTHfYTH2fANS2) 

CALL SQUAD 1(2fNTHfTH1fTH2fXTHfYTH3fANS3) 

FM 1 =FM1+ANS 1 
FH2=FM2+ANS2 
FZ 1 =FZ 1 + ANS3 

IF ( I WRT . EQ ♦ 1 ) WRITE( 11 f 110) TH 1 f TH2 f ANSI f 
1 ANS2» ANS3fYTH1 f YTH2f YTH3 
40 CONTINUE 


C 


FM1C=.5*FKR0C*RSP*RSP*GAM 
FM 1=FM1 -FM 1 C 
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135: 

136: 

137:C 

138: 

139: 

140: 

i4i :c 

142: 

143 J 

1 44 : c 
145 : 
146: 
1475C 
1 48 : c 
1 49 : c 
iso: 

151 J 

152: 

153:c 

154: 

155 : 

156: 

157 : c 

158: 

159: 

i6o:c 

i6i : 
162: 
163 : 
164: 

165 : c 
166 : 
167: 
168: 
i69:c 
170: 

171 : 

172; 

i73:c 

i74:c 

175: 

176: 


FMX=FM1*CPSI-FM2*SPSI 

FMY=FM1*SPSI+FM2*CF'SI 

FC1=FKR0C*CZ 

FC2=BETA*<P0-PI )*Rlhl/(2.+BETA*RIMl ) 
FC3=P0+PI+FC2 

FC=FC1+PIE*RM*RIM1*FC3 
FZ=FZ 1 “FC 

IF ( I WRT * GE 4 1 ) WRITE<11f112) FMX f FH Y f FZ f FC f FM1 f 
1 FM2fPSIfPSIDfFZ1 fANS1fANS2fANS3 


102 FORMAT ( 1H0 f 7 FQD2 NR 7 I2f 7 RI= 7 E12.5/ 

1 4E14*5/(4E14*5) ) 

104 FORMAT ( 1H0 f 7 FQD4 J4 7 I2f 7 NTH 7 1 2 f 3E15 . 5/ 7 X 7 5F8 *4/ 7 Y 7 5F8,4> 

108 FORMAT (1H0» 7 FGD8 NTH 7 fI3f 7 THA= 7 E 1 2 ♦ 5 » 7 STHA= 7 f E 1 2 . 5 f 

1 7 CTHA= 7 f E12 ♦ 5 f 7 ANSI = 7 f E12 ♦ 5 f 7 ANS2= 7 f E 1 2 . 5 f 

2 7 ANS3= 7 E12*5f/» ( 4E 1 5 ♦ 5 ) ) 

110 FORMAT (IHOf 7 F R 1 0 7 f 7 TH1= 7 E12«5f 7 TH2= 7 E12«5r 

1 7 ANS1= 7 E12 * 5 f 7 ANS2= 7 E 1 2 ♦ 5 f 7 ANS3 = 7 E 1 2 . 5 f / ( 5E15 ♦ 5 ) ) 

112 FORMATUHOf 7 FMX 7 f E 1 2 ♦ 5 f 7 FMY 7 f E 12 * 5 f 7 FZ 7 fE12*5f 

1 7 FC 7 El 2 « 5 f 7 FM 1 7 f E 1 2 ♦ 5 f 7 FM2 7 fE12.5f 7 PSI 7 fE12.5f 

2 7 PSID 7 fE12*5f/f 7 FZ 1 7 E 1 2 # 5 f 7 ANS1 7 E12.5f 

3 7 ANS2 7 E 1 2 ♦ 5 f 7 ANS3 7 E12*5> 

106 FORMAT ( 7 FQD6'f 7 R= 7 E12*5f 7 R2= 7 fE12.5f 7 ECTHA= 7 E 1 2 ♦ 5 f 

1 7 H“ 7 E 1 2 ♦ 5 f 7 HI = 7 E12 * 5/ 7 H0= 7 E12*5r 7 HM= 7 E12*5 f 

2 7 A= 7 E12 # 5 f 7 PH= 7 E12*5f 7 PS= 7 E12*5» 7 P= 7 E12*5) 

116 FORMAT ( 1 HI 7 CAX 7 fE12*5f 7 CAY 7 E12,5f 7 CZ 7 fE12.5f 

1 7 CAXD 7 fE12*5f 7 CAYD 7 fE12#5f 7 CZD 7 fE12*5f 

2 7 GAM 7 f E 1 1 ♦ 4 f 7 GAMD 7 f E 1 1 » 4 ) - 


RETURN 

END 
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6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
2 2 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


S85C 

59:C 

60S 


SUBROUTINE SQUAD1 (HOD 

NUMERICAL INTEGRATION BY 

DIMENSION A ( 70 ) , H < 70 ) 
DOUBLE PRECISION A , H 
DATA <A(I),H<I)»I=1,2 
1 7. 7 459666924 1483D-01 
1 8. 888888888 88889D-01 
1 3. 39981043584856 D -01 
1 2. 369268850561 89 D -01 
1-0. D-00 

1 1. 7 1 324 4923791 7 OD -01 
1 2. 3861 91 860831 9 7Ii -01 
1 1 .29484966168870D-01 
1 4. 05845151 377397D-01 
1 4. 17959183673469 D -01 
1 7. 96666477 41 3627D-01 
1 3. 137 06645877887 D -01 
1 9 . 68160239507626 D -01 
1 1 .80648160694857D-01 
1 3. 24253423403809 D -01 
1 3. 30239355001260 D -01 
1 8. 65063366688985 D -01 
1 2,19086362515982D-01 
1 1 .48874338981631D-01 
DATA ( A ( I ) » H ( I ) » 1=29 , ! 
1 9. 78228658146057 D -01 
1 1. 25580369464905 D -01 
1 5. 19096129206812 D -01 
1 2. 62804544510247 D -01 
1 9. 81560634246719 D -01 
1 1 .06939325995318D-01 
1 5. 87317954286617 D -01 
1 2.33492536538355D-01 
1 9.84183054718588D-01 
1 9. 21214998377280 D -02 
1 6. 423493394 40340 D -01 
1 2. 07816047536889 D -01 
1-0. D-00 

1 3.51194603317520D-02 
1 8. 272013 15069765D-01 
1 1, 57203 1 67 1 581 94 D -01 
1 3. 19112368927890D-01 
1 2. 15263853463158D-01 
1 9. 37273392400706 D -01 
DATA (A(I)fH(I), 1=57 , 

1 8.48206583410427D-01 
1 1 .39570677926154D-01 
1 3, 94151347077563 D -01 
1 1.984314853271 12 D -01 
1 9.89400934991650D-01 
1 6.22535239386480D-02 
1 7.55404408355003D-01 
1 1 .49595988816577D-01 
1 2.81603550779259D-01i 
1 1 .894506104550680-01/ 


»N,XI»XO»X,Y»ANSUR) 
GAUSS FORMULA 
X< N ) » Y ( N ) 

)/ 

5.555555555555560-01, 
8. 61136311594053 D -01 , 
6.52145154862546D-01 , 
5, 3846931 01 05683 D -01 , 
5.68888888888889D-01 f 
6.61209386466265B-01, 
4. 6791393457 26 91 D -01 , 
7,41531185599394 D-OIf 
3. 81 830050505 1 1 9D-0 1 f 
9. 60289856497536 D-01 f 
2. 2238 10344533 7 4D-01 f 
1 ,83434642495650 D-01 f 
8 . 12743883615740D-02f 
6. 133714 32 70059 0D-01 f 
3.123470770400030-01, 
9. 7390652851 7172 D -01 f 
1.49451349150581D-01, 
4 , 333953941292470 - 01 , 
2. 95524224 7 14753D-01/ 
>)/ 

5.56685671 161 740D-02, 
7. 30 15200557 4049 D-01 f 
2. 331 9376459 1990D-01 f 
0. D-OOf 

4. 7 17533638651 20D-02f 
7. 69902674 194305D-01f 
2. 03167426723066 D-01, 
1 .2523340851 1469D-01, 
4. 04840047653160 D -02, 
8.01578090733310D-01, 
1 .78145980761946D-01 , 
2. 3045831 59551 35D-01 , 
2. 32551553230874 D-01 , 
9.284348B3663574D-01, 
1.21518570687903 D-01, 
5. 152486363581 54 D-01, 
2. 05 19846372 1296D-01, 
9. 8799251 8020485D -01 , 
7. 0366047488 1080D-02/ 
)/ 

1 .07159220467172D-01 , 
5.70972172608539D-01, 
1.86161000015562D-01, 
0. D-00, 

2.71524594117540D-02, 
8.65631202387832D-01, 
1 .2462897 1255534D-01, 
4. 58016777657227 D-01, 
1.82603415044924D-01, 


S0U00020 

SGU00030 

SQU00040 

SGU00050 

SGU00060 

-0. D-00,SGU00070 

3. 4 78548451 374540-0 1 , SGU00080 
9. 06 179845938664 D-01fSGU00090 
4. 786286704993660 -01»SQU00100 
9. 324695142031 52D-01 fSGUOOI 10 
3. 6076 1 5730481 390-01 , SGU00120 
9. 4910791 2342759D-01 fSGUOO 130 
2. 797053914892770 -01, SGU00140 
-0. D-00 , SGU00150 

1 .012285362903760-01 ,SGUOO 160 
5. 255324099 1 6329 D-01,SQU00 170 
3.62683783378362D-01 ,SGU00180 
8. 3603 11 07 326636 D-01 ,S0U00190 
2.60610696402935D-01 , SQU00200 
-0. 0-00,SGU00210 

6. 6671 3443086880 D -02, SDU00220 
6. 794095682990240 -01, SGU00230 
2. 6926671930999 6D-01 ,S0U00240 

SQU00250 

SGU00260 

8. 87 062599768095 D-01,SQU 00270 
1 . 862902 10927734 D-01, SGU00280 
2. 6954 3 155952345 D-01 ,SQU00290 
2. 7292508677790 1D-01,SGU00300 
9.04U7256370475D-01 ,SGU00310 
1 ,60078328543346D-01,SGU00320 
3. 67831498998180 D -01, SGU00330 
2.49147045813403D-01 ,SGU00340 
9. 17598399222978D-01,SQU00350 
1 . 3887351 021 9787D-0 1 , SGU00360 
4. 48492751036447 D -01, SGU00370 
2.26283180262897D-01 ,SGU00380 
9. 862838086968 1 2D-01 , S0U00390 
8.01580871597600D-02»SQU00400 
6. 872929048 1 1685D-0 1 ,SGU004 10 
1 ,85538397477938D-01,SGU00420 
1.08054948707344D-01 ,SGU00430 
3. 07532419961 1 70D-02 , SGU00440 

SQU00450 

SQU00460 

7. 2441 7731360 170D-01 , SQU00470 
1 . 66269205816994D-01 ,SGU00480 
2.01194093997435D-01 ,SGU00490 
2. 0257824 1 925561D-01 ,SGU00500 
9.44575023073233D-01 ,SQU00510 
9.51585116824930D-02,SGU00520 
6. 17876244402644D-01 ,S0U00530 
1 .69156519395003D-01 ,SGU00540 
9.50125098376370D-02,SQU00550 

SGU00560 


NN = N 


61 

62 

63 

64 

65 

66 
67 


X IN= XI 
ANS=0 . 

IF ( X IN . EG . XO > GO TO ( 400 , 500 ), MODE 
IF (NN.LE. 16) GO TO 202 
WRITE (6,200) NN 

200 F0RMAT<47H0*****NUMBER OF POINTS FOR SGUAD1 INTEGRATION =I4»22H IS 
A OUT OF RANGE.*****) 
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68: 


RETURN 

69 J 

202 

INDKT = MOD ( NN j 2 > + 1 

70SC0MhENT--INDKT=l » N IS EVEN 

7i :c 


=2> N IS ODD 

72 : 


GO TO (204*210) t INDKT 

73 : 

204 

M IN = (NN*NN)/4 - 1 

74 : 


MAX = <NN*(NN+2> )/4 - 2 

75 : 


GO TO 215 

76: 

210 

M I N= <NN*NN-9)/4 + 1 

77 : 


MAX = (NN*(NN+2)-ll)/4 

78: 

215 

MINM1 = MIN-1 

79 J 


HDELX= 0 ♦ 5* ( XO-XIN ) 

so: 


M = NN/2 

8 1 : 


DO 220 1 = 1 f M 

82 : 


IM= I + M 

83: 


J= MINM1+I 

84: 


GO TO ( 216 r 217 > r MODE 

85: 

216 

X(I)= XIN + HDELX* < 1 *-A ( J ) ) 

86: 


X ( I M ) = XIN + HDELX* ( 1 . + A < J ) ) 

87: 


GO TO 220 

88: 

217 

ANS= ANS + (Y(I)+Y(IM) )*H< J> 

89: 

220 

CONTINUE 

90 : 


IF ( I NDKT ♦ EG * 2) GO TO ( 230 * 240) * MODE 

9i: 


GO TO (501 r 500) r MODE 

92 : 

230 

X ( NN ) = XIN+HDELX 

93 : 


RETURN 

94 : 

240 

ANS = ANS + Y(NN)*H(MAX+1 ) 

95 : 

500 

ANSUR= HDELX*ANS 

96: 

501 

RETURN 

97 : 

400 

DO 402 1 = 1* NN 

98: 

402 

X(I)= XIN 

99 : 

999 

RETURN 

loo: 


END 
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it 

2 ♦ C 

3JC 

4:c 

s:c 

6:c 

7: 

QIC 

9: 

io:c 

11: 

12: 

i3:c 

1 4 : 

15: 

i6t 

17: 

I8:c 

i9:c 
20: 
21 : 

22IC 

23: 
24: 
25: 
26: 
27: 
28: 
29: 
30: 
31 : 
32: 

33JC 

34: 

35: 

36 : c 

37 : c 

38:c 

39: 

40: 

4i:c 

42: 

43: 

44! 

45: 
46 J 
47:c 
48: 

49: 
so:c 
51 : 

52! 
53 :c 
54:c 
55 : 

56: 

57: 

58 : c 
59:c 
60: 
6i : 

62 :c 

63: 

64: 


SUBROUTINE WRT < N6 f JT f Y f T f DELH f TEND f INDEX f IER > 


THIS SUB WRITES OUTPUT FOR ETSION'S PROGRAM 


1 


DIMENSION Y(N6) 

COMMON/COM1 / PIE f PIE2 f N4FIE f NR f NTH f IWRT 

C0MM0N/C0M2/ FKROC f RQ f BETA f RI »COf 
RMfPOfPIfFI fFMEfFMfRSP 


C0MM0N/C0M3/ HMIN f HMIN1 f HMIN2 
COMMON/COM4/ CAX f CAY f CZ f C AXD f CAYD f CZD 


COMMON/COM5/ FMX f FMY f FZ f GAM f PSI f GAMD f PSID f FC f 
1 FM1 fFM1CfFZ1fFM2fEPSfDEL 


TURT = TEND 

IF ( JT • LE • 1 ) TWRT = 0 « 


CAX= Y ( 1 ) 

CA Y = Y < 2 ) 

CZ = Y < 3 ) 

CAXD=Y (4) 

C A Y D = Y < 5 ) 

CZD=Y<6> 

CALL FGUAD 
YPRSV1=(FMX+FME)/FI 
YPRS V2=FM Y/F I 
YPRSU3=FZ/FM 

WRITE (7 f 102 ) INDEXfIERfY( 4 ) f Y ( 5 ) f Y <6 ) f T f DELH f YPRSV 1 f YPRSV 2 f YPRSV 3 f 
1 FC 


SGXY=CAX*CAX + CAY*CAY 
GAM= SORT ( SQXY ) 

CAXX = ABS(CAX) 

CAYY= ABS(CAY) 

PSI= ATAN2(CAYYfCAXX) 

IF ( CAX.LT ♦O.AND.CAY.GT*0) PSI=-PSI 
IF(CAX.GE.O.AND.CAY.LT.O) PSI=-PSI 

CPSI =CAX/GAM 
SPS I =C AY/GAM 

GAMD=CAXD*CPSI+CAYD*SPSI 

PSID=<CAYD*CPSI-CAXD*SPSI)/GAM 


TWRTP=TWRT/PIE2 

WRITE(6 f 104 ) CZ f GAM f PSI f CZD f GAMD f PSID t CAX f 
1 CAYfHMINfTWRTP 


104 FORMAT < 1H rlOE11.4> 
102 FORMAT (2I5 f9E11 *4) 

90 RETURN 
END 
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APPENDIX D - FLOW CHARTS 


MAIN EXECUTIVE PROGRAM 
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37 
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APPENDIX E 
SAMPLE PROBLEM 


A computer output listing is shown as an example of a representative run. 
The computer run generated output for twenty revolutions but only output from 
the first, second, nineteenth, and the twentieth revolutions is shown. Input 
data of the two namelists, INPUT1 and INPUT2 is shown preceeding the listed 
output and this input is defined in the report section, USING THE PROGRAM. 


In this case, the Adams predictor corrector method was used (METH = 1) , an 
initial step size of 0.001 was tried (DELH = 0.001), functional iteration was 
used (MITER = 0), and the input relative error bound, T0L, was set to 0.0001. 

The initial values T = 0 (TZR0 = 0) and a x = 0, a y = 0.1, Z = 0, a x = 0, 

oty = 0, and Z = 0 (YIN = 0, 0.1, 0,0, 0,0) were set. TSTOP was set equal to 

125.66 radians, or twenty revolutions of the rotor. Also, the inputs R^ = 0.9, 

C o 

K — = 0.0025, B = 10, P ± = 0, P Q = 0.25, Mg = 0, m = 0.08, r g = 1, and RSP = 1 
were used. 


The execution of this case required two minutes of computer time on the 


Univac 

1100/42. 

SlNPUTl 
READ 1 


F 


METH 

z 


♦1 

MITER 

r 


♦ 0 

T0L 

r 


.10000000E-Q3 

INDEX 

z 


♦1 

DELH 

— 


•10000000E-02 

T2R0 

z 


•OOOOOOOOE^OO 

YIN 

r 


•OOCOOOOOE*QO 

TSTOP 

r 


.OOOOOOOOE^OO 
• 12566 3 71 E *03 

DTEN0 

z 


.1 2 566 3 71 E ♦00 

NR 

z 


♦ 5 

NTH 

z 


♦ 8 

HMST0P 

z 


• 100000 00E *20 

HS TAP 

z 


• 1 0000000E -02 

NPLT 

z 


♦1 

PUNCH 

z 

F 


SEND 
SINPUT2 
R 1 



• 9 OOOOOOOE -*00 

FK ROC 

z 


•250000 00 E ~02 

BETA 

z 


• 1 00000 00 E +02 

PI 

z 


•OOOOQOOOE*Oo 

P0 

z 


•2SOOOOOOE+QO 

FME 

z 


•OOOOOOQOE 4 QO 

FM 

z 


• 800000 OOE— 01 

RG 

- 


•1000D000E+Q1 

RSP 

- 


• 1 OOOOOOOE 


. 10 0000 0 QE 4 00 » • 000 000G0E *00 , • 000Q0000E *00 , 

• 00000000E»00 


SEND 



cz 

BAH 

PSI 

.0000 

. 1000*00 

. 1571*01 

. 1153-03 

. 9986-01 

-. 1569*01 

.* 605-03 

• 99**-0 1 

-. 1565*01 

• 1021 -02 

. 9876-01 

-. 1559*0 1 

. 178* -02 

. 978* -0 1 

-. 1550*0 1 

. 2735-02 

. 9670-01 

- • 1538 <*0 1 

.3871 -02 

. 9535-01 

-. 152**0 1 

. 5150-02 

. 9383-01 

-. 1507*0 1 

. 6590-02 

*92 1*“0 1 

-. 1 * 87*0 1 

.8171 -02 

. 9030-01 

-. 1 * 65*0 1 

. 9878-02 

. 8835-01 

1 ** 0«0 1 

. 1169-01 

. 8631-01 

-. 1 * 13*0 1 

. 1361-01 

. 8 * 20-01 

-. 1382*0 1 

. 1560-01 

• 820*-0 1 

-. 13 * 9*0 1 

. 1768-01 

. 7985-01 

-. 1312*01 

. 1982 -Ql 

. 7766-01 

-. 1273*0 1 

. 220* -01 

. 75 * 9-01 

-. 1230*0 1 

• 2 * 26-01 

. 7336-01 

-. 1183*0 1 

. 265* -01 

. 7129-01 

-. 113**01 

. 2885-01 

. 6930-01 

-• 1 Qfl 1*0 1 

. 3117-01 

. 67 * 2-01 

-• 102 * *0 1 

. 3351-01 

• 6567-01 

-. 96 * 8*00 

. 3585-01 

• 6 * 06-01 

-. 9019*0 C 

. 3819-01 

. 6261-01 

-. 8359*00 

.* 051 -01 

. 6133-01 

-. 7673*00 

.* 281-01 

. 602 *- 0 l 

-. 6962*00 

.* 509-01 

. 5936-01 

-. 6231*00 

.* 73* -01 

. 5867-01 

-. 5 * 85*00 

.* 956-01 

. 5819-01 

-.* 730*00 

. 5173-01 

. 5791-01 

-. 3969*00 

. 5386-01 

. 5783-01 

-. 3209*00 

. 559* -01 

. 5793-01 

-. 2 * 55*0 0 

. 5797-01 

. 5820-01 

-. 1710*0 0 

. 5993-01 

. 5862-01 

- . 9 7 8 * -0 1 

. 6183-01 

. 5917-01 

-. 2627-0 1 

• 6367-01 

. 5985-01 

.*35 0-0 1 

. 65 * 3-01 

. 6063-01 

. 111**0 c 

. 6713-01 

. 6150-01 

. 1772*00 

• 6876-01 

. 62 * 2-01 

. 2 * 11*00 

. 7030-01 

. 6339-01 

• 30 3 1 *Q C 

. 7177-01 

• 6 * * 0-01 

.36 33*00 

. 7317-01 

. 65 * 3-01 

• * 2 17*0 0 

. 7 * 1 * 8-01 

• 66 * 7-01 

.* 785*00 

. 7572-01 

. 6751-01 

. 5338*00 

. 7687-01 

. 6853-01 

. 5877*00 

• 779* -01 

. 6953-01 

• 6 * 05*00 

. 7893-01 

. 7050-01 

• 6922*00 

• 798* -01 

. 71 * 3-01 

. 7 * 29*00 

. 8067-01 

. 7233-01 

. 7928*00 

. 81 * 2-01 

. 7318-01 

. 8 * 20*00 

. 8208-01 

. 7399-01 

. 8907*0 0 

. 8267-01 

• 7 * 7*-0 1 

. 9389*00 

. 8318-01 

• 75**-0 1 

. 9867*0 0 

. 8362-01 

. 7609-01 

• 103**0 1 

. 8399-01 

. 7669-01 

• 1082*0 1 

• 6 * 28-01 

. 772*-0 1 

. 1129*0 1 


CZO 

G A PD 

PS 10 

.0000 

.0000 

• 0000 

. 1881-02 

-. 2279-02 

. 2177-01 

. 36 5* -02 

-.** 22-02 

.* 282-01 

. 5318-02 

-. 6 * 21-02 

. 6338-01 

. 6873-02 

-. 8268-02 

• 8 367 - Cl 

. 8320-02 

-.9957 -02 

• 1038*00 

. 9658-02 

-. 11 * 8-01 

• 12 * 1 *00 

. 1089-01 

-• 1 28* -01 

• 1 * * 5 ♦ 00 

. 1202-01 

-. 1 * 03-01 

. 1653*00 

. 13 G 5-01 

-. 150* -01 

• 1 866 ♦ 00 

. 1399-01 

-. 1587 -01 

. 2085*00 

. 1 * 82-01 

-. 1653-01 

. 2310*00 

. 1556-01 

-. 1 701 -01 

. 25 * 1*00 

. 1622-01 

-.1 732-01 

. 2780*00 

. 1678-01 

-. 17 * 5-01 

. 3027*00 

. 1727-01 

-. 1739-01 

. 3282*00 

. 1767-01 

-. 1716-01 

. 35 * 3*00 

. 1800-01 

-.1675 -01 

. 3810*00 

. 1825-01 

-. 1616-01 

• 4080*00 

. 18 * 3-01 

-. 15 * 1-01 

•* 351*00 

. 1855-01 

—. 1 ** 9 -01 

.* 619*00 

. 1860-01 

-. 13*1 -01 

.* 880*00 

. 1860-01 

-. 1219-01 

. 5128*00 

. 1853-01 

-. 1 QS* -Ql 

. 5358*00 

. 18 * 1-01 

-. 9382-02 

.55 6 ** CO 

• 182* -01 

-. 78 * 5-02 

. 5739*00 

. 1802-01 

-. 625* -02 

. 5878*00 

• 1776-01 

-.* 6 * 1-02 

. 5978*00 

. 17 * 5-01 

-. 3037-02 

. 6037*00 

. 1710-01 

-. 1 * 73-02 

.6053 *00 

• 1671-01 

.2301 - 0 * 

. 6029*00 

. 1629-01 

. 1 * 2* -02 

. 5969*00 

• 15 8 *-Ql 

. 2710-02 

. 5879*00 

• 1536-01 

. 3866 -02 

. 5763*00 

• 1 * 85-01 

.* 882-02 

. 5629*00 

. 1 * 32-01 

. 5755-02 

. 5 * 82*00 

. 1376-01 

. 6 * 85-02 

. 5327*00 

. 1319-01 

. 7074-02 

. 5170*00 

. 1260-01 

. 7530-02 

. 5015*00 

. 1199-01 

. 7863-02 

.* 865*00 

. 1137-01 

• 8080-02 

.* 723*00 

. 107*-01 

. 8192-02 

.* 589*00 

. 1011-01 

. 8209 -02 

. 4*65 *00 

. 9 * 67-02 

. 81 ** -02 

. 4352*00 

. 8822-02 

. 8006-02 

. 4251*00 

• 6176-02 

. 7806-02 

. 4160*00 

. 7530-02 

. 7553-02 

. 4080*00 

. 8885-02 

. 7254-02 

. 4010*00 

. 62 * 5-02 

. 6918-02 

. 3950*00 

. 5610-02 

. 6553-02 

. 3899*00 

.* 983-02 

.6 166 -02 

. 3857*00 

.* 365-02 

. 576* -02 

. 3823*00 

. 3757-02 

. 5351-02 

. 3797*00 

. 3161-02 

. 4934-02 

. 3778*00 

. 2578-02 

, 4516-02 

. 3766*00 

. 2010-02 

. 4 103-02 

. 3759*00 


C A X 

cax 

HHTN 

.0000 

. 1000*00 

.913 G *00 

-. 1335-03 

. 9986-01 

.0 101*00 

-. 5332-03 

. 99 * 4-01 

. 9105*00 

-. 1183 - 0 ? 

. 9876-01 

. 9112*00 

-. 2072-02 

. 9782-01 

• 9121*00 

-. 3180-02 

. 9665-01 

• 9132*00 

-.* 508-02 

. 952**01 

. 91 * 5*00 

-. 5999-02 

. 9363-01 

. 9160*00 

-. 7676-02 

. 9182-01 

. 9176*00 

-. 9516-02 

. 8980-01 

. 9194*00 

-. 11 * 9-03 

. 8760-01 

. 9213*00 

-. 1359-01 

. 8523-01 

. 9232*00 

-. 1579*01 

. 8271-01 

. 9252*00 

-. 1807-01 

• 8003-01 

. 9273*00 

-. 20 *» 2 - 0 l 

. 7720-01 

. 929*400 

-. 2202-01 

. 7 * 23-01 

. 9315*00 

-. 2525-01 

.7 1 1*-01 

. 9335*00 

-. 2771-01 

. 6792-01 

. 9355*00 

-. 3017-01 

. 6 * 59-01 

. 9375*00 

-. 3261-01 

. 6115-01 

. 939**00 

-. 3503-01 

. 5761-01 

• 9 * 12*00 

-. 37 * 1-01 

. 5398-01 

. 9 * 28*00 

-. 3973-01 

. 5025-01 

. 94 * 3*00 

-.* 198-01 

.* 6 * 5-01 

. 9 * 57*00 

-.** 15-01 

.* 257-01 

. 9470*00 

-. 4623-01 

. 3863-01 

• 9480*00 

-.* 820-01 

. 3 * 64-01 

. 9489*00 

-. SQQ 6-01 

. 3059-01 

• 9496*00 

-. 5180-01 

. 2651-01 

. 9501*00 

-. 53 * 1-01 

. 2239-01 

. 950**00 

-. 5 * 87-01 

. 1824-01 

• 9506*00 

-. 5619-01 

. 1 * 08-01 

• 9506*00 

-. 5735-01 

. 9904-02 

. 9505*00 

-. 5834-01 

. 5726-02 

. 9502*00 

-. 5915-01 

. 1555-02 

. 9498*00 

-. 5979-01 

-. 2603-02 

. 949**00 

-. 6026-01 

-. 6737-02 

. 9488*00 

-. 6053-01 

-. 1084-01 

. 9 * 81*00 

-. 6062-01 

-. 1 * 91-01 

. 947**00 

-. 6050-01 

-. 1892-01 

. 9 * 67*00 

-. 6020-01 

-. 2289-01 

. 9459*00 

-. 5970-01 

-. 2678-01 

. 9 * 51*00 

-. 5901-01 

-. 3061-01 

. 9 * 43*00 

-. 5812-01 

-. 3 * 35-01 

. 9435*00 

-.5 703-01 

-. 3800-01 

. 9 * 27*00 

-. 5575-01 

-. 4155-01 

. 9 * 20*00 

-. 5 * 27-01 

-. 4 * 99-01 

. 9 * 12*00 

-. 5261-01 

-. 4832-01 

. 9 * 05*00 

-. 5076-01 

-. 5152-01 

. 9398*00 

-. 4874-01 

-. 5 * 60-01 

. 939 1*00 

-.465 3-01 

-• 5 7 S 3 -01 

. 9385*00 

-. 4 * 15-01 

-. 6031-01 

. 9379*00 

-. 4160-01 

-. 6293-01 

. 9373*00 

-. 3890-01 

-. 65 * 0-01 

. 9368*00 

-. 3604-01 

-. 6770-01 

. 9363*00 

JJ 05-01 

-. 6982-01 

. 9359*00 


T MR T 
• 0000 
. 2000-01 

• «|QOO-Ol 
. 6000-01 
. 8000-01 
. 1000*00 
. 1200*00 
* 1 * 400*00 
. 1600*00 
. 1600*00 
. 2000*00 
. 2200*00 
. 2 * 400 * 00 
. 2600*00 
• 2800*00 
. 3000*00 
. 3200*00 
. 3*4 00*00 
. 3600*00 
• 3800*00 
.* 000*00 
. 41200*00 

• ** 400*00 
.* 600*00 
. 41800*00 
*500 o * 00 
. 5200*00 
. 5 * 00*00 

. 5600*00 O 

. 5800*00 
• 6000*00 
. 6200*00 

• 6 * 400*00 
•6 Cqq* 00 
. 6800*00 
. 7000*00 
. 7200*00 
. 7 * 00*00 
. 7600*00 
. 7800*00 
• 8000*00 

• 8200*00 
. 8 * 00*00 
. 8600*00 
. 8800*00 
. 9000*00 
. 9200*00 
• 9 * 00*00 
. 9600*00 
. 9800 * 00 
• 1000*01 
.102 0*01 
. 10 * 0*01 
. 1060*01 
• 108 0*01 
. 1100*01 



. 8450-01 

. 7774-01 

» 11 76 * 0 1 

. 1458-02 

.3699 -02 

.3758 *00 

-. 2991-01 

-. 7375-01 

. 9355*00 

. 1120*01 

. 8465 -01 

. 7818-01 

• 122 3*0 1 

. 9240-03 

.3 306-02 

. 3762*00 

-. 2664-01 

-. 7350-01 

. 9351*00 

.1 140*01 

. 8473-01 

. 7857-01 

. 1270*0 1 

. 4074-03 

.2929 *0 2 

. 3771*00 

-. 2325-01 

-. 7505-01 

. 9348*00 

. 1160*01 

. 8475-01 

. 7892-01 

. 1318*01 

-. 9053-04 

.2569 -02 

. 3783*00 

-. 1976-01 

-. 7641-01 

. 9345*00 

. 1180*01 

.8471 -01 

. 7923-01 

• 1365*01 

-. 5687-03 

. 2230-02 

. 3800*00 

-. 1617 -Ui 

-. 7757-01 

• 5 34 3 * QQ 

. 1200*01 

.8461 -01 

. 7950-01 

• 1413*01 

-. 1026*02 

. 1913-02 

. 3820*00 

-. 1248-01 

-. 7851-01 

. 9340*00 

. 1220*01 

. 8446-01 

. 7972-01 

• 1461*01 

-. 1462-02 

. 1620-02 

. 3842*00 

-. 8720-02 

-. 7925-01 

. 9338*00 

* 124 { j * ol 

. 8425-01 

. 7991-01 

« 1 5 1 Q*D 1 

-. 1874-02 

. 1352-02 

. 3867*00 

-. 4887-02 

-. 7976-01 

. 9337*00 

. 1260*01 

. 8399-01 

. 8007-01 

. 1558*0 1 

-. 2265-02 

. 1111-02 

. 3894*00 

-. 9984-03 

-. 0007-01 

. 9335*00 

• 1280*01 

. 8368-01 

• 8020-01 

-. 1534*01 

-. 2633-02 

. 8979-03 

. 3923*00 

. 2934-02 

-. 8015-01 

. 9334*00 

• 1 3 Q 0 * 01 

. 8333-01 

. 8031-01 

-. 1485*0 1 

-. 2978-02 

. 7121-03 

. 3953*00 

. 6898-02 

-. 8002-01 

. 9333*00 

. 1320*01 

.8294 -01 

. 8040-01 

-. 1435*01 

-. 3298-02 

. 5543*03 

. 3984*00 

. 1008-01 

-. 7966-01 

. 9332*00 

• 1340*01 

. 8250-01 

. 8046-01 

-. 1385*0 1 

-. 3595-02 

. 4242-03 

. 4015*00 

. 1488-01 

-. 7908-01 

. 9331*00 

. 1360*01 

. 8204-01 

. 8051-01 

-. 1334*01 

-. 3866-02 

. 3213-03 

. 4047*00 

. 1887-01 

-. 7827-01 

. 9330*00 

• 1380*01 

. 8155-01 

. 8055-01 

-. 1283*0 1 

-. 4114-02 

. 2453-03 

. 4078*00 

. 2285-01 

.-. 7724-01 

• 9330*00 

• 1400*01 

. 8100-01 

. 8058-01 

-. 1232*0 1 

-. 4338-02 

. 1954-03 

• 4110*00 

. 2680-01 

-. 7599-01 

. 9329*00 

. 1420*01 

• 8044-01 

. 8061-01 

-. 1180*01 

-. 4538-02 

. 1702*03 

• 4140*00 

. 3071-01 

-. 7453-01 

. 9329*00 

. 1440*01 

. 7986-01 

. 0063-01 

-. 1128*01 

-. 4714-02 

• 1664 *03 

. 4170*00 

. 3457-01 

-. 7284-01 

• 9328*00 

. 1460*01 

. 7926-01 

• 8066-01 

-. 1075*01 

-. 4867-02 

.1 885 *03 

. 4199*00 

. 3836-01 

-. 7095-01 

. 9327*00 

. 1480*01 

. 7864-01 

. 8069-01 

-. 1022*01 

-. 4997-02 

. 2292-03 

. 4226*00 

. 4208-01 

-. 6885-01 

. 9327*00 

. 1500*01 

. 7800-01 

. 8072-01 

-. 9690*00 

-. 5104-02 

. 2886-03 

. 4251*00 

. 4570-01 

-. 6654-01 

. 9326*00 

. 1520*01 

. 7736-01 

. 8077-01 

-. 9154*00 

-. 5189-02 

.3647 -03 

. 4275*00 

. 4922-01 

-. 6403-01 

. 9325*00 

• 1540*01 

. 7670-01 

. 8082-01 

-. 8616*00 

-. 5252-02 

. 4552-03 

. 4296*00 

. 5263-01 

-. 6133-01 

• 9324*00 

. 1560*01 

. 7604-01 

. 8089-01 

-. 8075*00 

-. 5294-02 

. 5585-03 

• 4316*00 

. 5592-01 

-. 5844—01 

. 9323*00 

• 1580*01 

. 7537-01 

. 8097-01 

-. 7532*00 

-. 5316-02 

.6725 -03 

. 4333*00 

. 5907-01 

-. 5538-01 

. 9322*00 

• 1600*01 

. 7470-01 

. 8106-01 

-. 6986*00 

-. 5319-02 

. 7951-03 

. 4349*00 

. 6207-01 

-. 5214-01 

• 932 1 * Q 0 

. 1620*01 

. 7405-01 

• 8118-0 1 

-. 6439*00 

-. 5303-02 

. 9240*03 

. 4362*00 

. 6492-01 

-. 4873-01 

. 9320*00 

• 1640*01 

.7337 -01 

. 8130-0 1 

-. 5890*0 0 

-» 52 b 8 - Q 2 

. 1057-82 

. 4372*00 

. 6760-01 

-. 4517-01 

. 9318*00 

. 1660*01 

. 7271-01 

• 8145-01 

-. 5341*00 

-. 5216-02 

. 1193-02 

. 4381*00 

. 7010-01 

-. 4146-01 

. 9317*00 

• 168 0*01 

. 7205-01 

. 8161-0 1 

-. 4790*00 

-. 5148-02 

. 1330-02 

. 4387*00 

. 7243-01 

-. 3761-01 

. 9315*00 

. 1700*01 

. 7141-01 

. 8179-01 

-. 4238*00 

-. 5065-02 

. 1466-02 

. 4392*00 

. 7455-01 

-. 3364-01 

. 9313*00 

. 1720*01 

. 7078-01 

. 8199-01 

-. 3687*00 

-. 4967-02 

. 1599-02 

. 4395*00 

. 7648-01 

-. 2954-01 

. 9311*00 

. 1740*01 

. 7016-01 

. 6220-01 

-. 3135*00 

-. 4855-02 

. 1729-02 

. 4395*00 

. 7820-01 

-. 2535-01 

. 9309*00 

. 1760*01 

. 6956-01 

. 8243-01 

-. 2582*00 

-. 4731-02 

. 1853*02 

. 4395*00 

. 7969-01 

-. 2105-01 

. 9306*00 

. 1780*01 

. 6897-01 

. 8267-01 

-. 2030*00 

-. 4594-02 

. 1972-02 

. 4392*00 

. 8097-01 

-. 1667-01 

. 9304*00 

• 1800*01 

. 6840-01 

. 8293-01 

-. 1479*00 

-. 4448-02 

. 2083 -02 

. 4389*00 

. 8202-01 

-. 1222-01 

. 9301*00 

. 1820*01 

. 6785-01 

. 8320-01 

-. 9278-0 1 

-. 4291-02 

. 2187-02 

.4384 *00 

. 8284-01 

-. 7708-02 

. 9299*00 

. 1840*01 

. 6732-01 

• 8348-01 

-. 3772-01 

-. 4126-02 

. 2282 -02 

. 4379*00 

. 8342-01 

-. 3148-02 

. 9296*00 

• 1860*01 

. 6682-01 

. 8378-01 

. 1728-0 1 

-. 3952-02 

. 2369-02 

. 4372*00 

. 8376-01 

. 1448-02 

• 9293*00 

. 1880*01 

• 6635 - Q 1 

• 8408-01 

. 7220-01 

-. 3772-02 

. 2446 -02 

• 4365*00 

. 8386-01 

.606 £,-02 

. 9290*00 

. 1900*01 

. 6587-01 

. 8439-01 

. 12 70*0 C 

-. 3586-02 

. 2515 -02 

. 4358*00 

. 8371-01 

. 1069-01 

. 9287*00 

. 1920*01 

. 6543-01 

. 8471-01 

• 1818*0 0 

-. 3395-02 

. 2575-02 

. 4350*00 

. 8332-01 

. 1531-01 

• 9284*00 

. 1940*01 

• 6 $Qt -01 

. 8504-01 

. 2364*00 

-. 3199-02 

. 2626 -02 

. 4342*00 

. 8267-01 

. 1992-01 

• 9281*00 

. 1960*01 

. 6463 -01 

. 8537-01 

• 2909*00 

-. 3001-02 

.2669 - Q 2 

. 4334*00 

. 8179-01 

. 2449-01 

. 9278*00 

. 1980*01 

. 6426-01 

. 8571-01 

. 3453*00 

-. 2800-02 

. 2704-02 

. 4325*00 

. 8065-01 

. 2901-01 

. 9275*00 

. 2000*01 



. 5320*00 

. 1490*01 

-. 8256*00 

. 5360*00 

. 1497*01 

-. 7720*00 

.5399 *00 

. 1504*01 

-. 7184*0 C 

. 5439*00 

. 1510*01 

-. 6649*00 

• 548D ♦ 00 

. 1517*01 

-.61 13*00 

. 5521 *00 

. 1524*01 

-. 5577*00 

.5562 *00 

. 1531*01 

-. 504 1*00 

.5603 *00 

. 1538*01 

-. 4505*00 

,5695 *00 

. 1545*01 

-. 3969*0 0 

.5688 *00 

. 1552*01 

-. 3433*0 C 

. 5731 *00 

. 1559*01 

-. 2897*00 

. 5774*00 

. 1566*01 

-. 2361*00 

• 5818*00 

. 1573*01 

-. 1825*0 0 

.5862 *00 

. 1580*01 

-. 1289*00 

.5907 *00 

. 1587*01 

-. 7535-01 

. 5952*00 

. 1594*01 

-. 2176-0 1 

. 5998*00 

. 1602*01 

. 3182*0 1 

.6014 *00 

. 1609*01 

. 8540-0 1 

.6090 *00 

.16 16*01 

. 1390*00 

. 6137*00 

.16 24*01 

. 1926*0 C 

. 6185*00 

• 1631*01 

• 2463*00 

. 6233*00 

. 1639*01 

. 2999*00 

. 6281*00 

• 1646*01 

. 3535*0 C 

• 6330*00 

• 1654*01 

. 4071*0 C 

. 6379*00 

• 1662*01 

. 4607*00 

.6429 *00 

• 1670*01 

. 5143*00 

. 6483 *00 

. 1677*01 

. 5679*00 

• 6531 *00 

• 1685*01 

• 6216*00 

• 6582*00 

. 1693*01 

• 6752*00 

.6634 *00 

. 1701*01 

. 7288*0 0 

. 6687*00 

. 1709*01 

. 7824*00 

. 6740*00 

. 1717*01 

* 8 36 Q«0 0 

.6794 *00 

. 1726*01 

. 8897*00 

. 6849*00 

. 1734 * 0-1 

. 9433*00 

.6904 *00 

. 17 * 42*01 

. 9969*00 

. 6960*00 

. 1750*01 

. 1051*0 1 

. 7016*00 

. 1759*01 

• 1104*01 

. 7072*00 

. 1767*01 

. 1158*01 

. 7130*00 

. 1776*01 

• 1211*0 1 

.7188 *00 

. 1784*01 

. 1265 * 0 1 

. 7247*00 

. 1793*01 

. 1319*0 1 

. 7306*00 

• 1802*01 

. 1372*0 1 

. 7366 * 0D 

. 1810*01 

. 1426*0 1 

. 7426*00 

. 1 819*01 

. 1479*0 1 

. 7488*00 

. 1826*01 

. 1533*01 

• 7550*00 

.18 37*01 

1555*0 1 

. 7612*00 

. 1846*01 

-. 1501*0 1 

. 7676*00 

. 1855*01 

-. 1448*0 1 

. 7740*00 

. 1865*01 

-. 1394*0 1 

. 7804 *00 

. 1074*01 

-. 1340*0 1 

. 7870*00 

. 1803*01 

-. 1287*01 

. 7936*00 

. 1893*01 

-. 1233*0 1 

. 8003*00 

. 1902*01 

-. 1180*0 1 

. 8070*00 

. 1912*01 

-. 1126*0 1 

. 8139*00 

. 1921*01 

-. 1072*0 1 

• 8208*00 

• 1931*01 

-. 1019*0 1 

• 8278 *00 

. 1941*01 

-. 96 5 1 *0 C 


. 3100-01 

.5 260-01 

. 4265*00 

. 3130-01 

. 5294-01 

. 4265*00 

. 3162-01 

. 5326-01 

. 4265*00 

. 3194-01 

. 5359-01 

. 4265*00 

. 3226-01 

. 5394-01 

. 4265*00 

. 3257-01 

.5431 -01 

. 4265*00 

. 3290-01 

. 5465-01 

. 4265*00 

. 3325-01 

. 5496 - U1 

. 4265*00 

. 3363-01 

. 5524 -01 

. 4265*00 

. 3397-01 

. 5560-01 

. 4265*00 

. 34 33-01 

• * 594 -01 

. 4265*00 

. 3465-01 

. 5634-01 

. 4265 ♦ 00 

. 3502-01 

• 5 669 -01 

. 4266*00 

• 3535-01 

. 5712-01 

• 4266*00 

. 3569-01 

. 5751-01 

.4 266 ♦ no 

• 3606-01 

. 5789-01 

.4266 *00 

. 3644-01 

. 5827-01 

• 4266 *CQ 

. 3682-01 

.5865 -01 

. 4266*00 

. 3720-01 

. 5905-01 

. 4266*00 

. 3758-01 

. 5946-01 

. 4266*00 

. 3796-01 

. 5908-01 

. 4266*00 

. 38 35-01 

. 6031 -01 

. 4266*00 

. 3874-01 

. 6073-01 

. 4267*00 

.39 15-01 

. 6114-01 

. 4267*00 

. 3958-01 

. 6153-01 

. 4267*00 

. 4002-01 

. 6192-01 

. 4267*00 

. 4047-01 

. 6230-01 

. 4267*00 

. 4093-01 

. 6269-01 

. 4267*00 

. 4139-01 

. 6309-01 

•4267400 

. 4184-01 

. 6352-01 

. 4267*00 

. 4228-01 

. 6399-01 

. 4267*00 

. 4271-01 

. 6449-01 

.4268400 

. 4313-01 

. 6502-01 

. 4268*00 

. 4356-01 

. 6557-01 

. 4268*00 

. 4402-01 

. 6607-01 

.4268 *QD 

. 4447-01 

. 6660-01 

. 4268*00 

. 4498-01 

. 6704 -01 

. 4268*00 

. 4548-01 

. 6753-01 

. 4268*00 

. 4598-01 

. 6801 -01 

. 4268*00 

• 4646-01 

. 6856-01 

• 4268*00 

. 4693-01 

• 6913-01 

. 4268*00 

. 4745-01 

. 6964-01 

. 4268*00 

. 4800-01 

. 7011-01 

. 4268*00 

. 4855-01 

. 7061-01 

. 4268*00 

. 4908-01 

. 7115-01 

• 4268*00 

. 4959-01 

. 7174-01 

. 4268*00 

. 5012-01 

. 7232-01 

. 4268*00 

. 5066-01 

.7291 -01 

.4268 *00 

. 5121-01 

. 7348 -01 

. 4268*00 

. 5180-01 

. 7402 -01 

. 4268*00 

. 5235-01 

. 7464 -01 

. 4268*00 

. 5290-01 

.7520 -01 

. 4268*00 

. 5348-01 

. 7588-01 

. 4268*00 

. 5409-01 

. 7645-01 

.4268 *00 

. 5472-01 

. 7702-01 

. 4268*00 

. 5536-01 

. 7758-01 

. 4268*00 

. 5601-01 

. 7815-01 

. 4268*00 


. 1011*01 
. 1073401 
.1132401 
.1 189*01 

• 1 29 2 * 01 
. 12*3401 
. 1340*01 
. 1384*01 
. 1925*01 

• 1461*01 

• 1 4*4401 
. 1522*01 
. 1547*01 
.1567401 
.1563401 
. 15 * 4*01 
. 1601*01 
. 1603*01 
. 1601*01 
. 15 * 4*01 
. 1582*01 
. 1566*01 
. 1545*01 
. 1519*01 

. 14 * 9 * 0 ! 

. 1454*01 
. 1414*01 
. 1370*01 
. 1322*01 
. 1269*01 
.1 212*01 
. 1151*01 
. 1087*01 
. 1018*01 
. 9458 ♦ 00 
.870 1*00 
. 7912*00 
. 7 0 * 3 4 00 
• 6245*00 
. 5372*00 
. 4474*00 
. 3555*00 
. 2616*00 
. 1661*00 
. 6911 - 0l 
-. 28 * 5-01 
-. 1280*00 
-. 2277*00 
-. 3277*00 
-. 4277*00 
-. 5275*00 
-. 6267*00 
-. 7252*00 
-. 8225*00 
-.9 184 * D0 
-. 1013*01 
-. 1105*01 


-. 1095*01 
-. 1044*01 
-. 9897*00 
-. 9318*00 
-. 8707*00 
-. 8064*00 
-. 7393*00 
-. 6695*00 
-. 5971*00 
-. 5222*00 
-. 4452*00 
-. 3662*00 
-.2854 *00 
-. 2031*00 
-. 1195*00 
-. 3469-01 
. 5095-01 
. 1372*00 
. 2239*00 
.3109400 
. 3977*00 
.484 1*00 
. 5700*00 
. 6550 * 00 
. 7388*00 
. 8214*00 
. 9323*00 
. 9814*00 
. 1058*01 
. 1133*01 
. 1205*01 
. 127^*01 
. 1341*01 
. 1404*01 
. 1463*01 
.1519401 
. 1571*01 
. 1619*01 
. 1662*01 
. 1702*01 
. 1736*01 
. 1766*01 
. 1791*01 
. 1812*01 
. 1827*01 
.1037401 
. 1842*01 
.184 1*01 
. 1836*01 
. 1824*01 
. 1808*01 
. 1786*01 
. 1758*01 
. 1726*01 
. 1687*01 
. 1644*01 
. 1595*01 
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